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ABSTRACT 



supernovae are 3.29+ 3 78 



"•^"-3.12-2.1 



'xlO -4 yr 1 



We present supernova rate measurements at redshift 0.1-1.0 from the Stockholm VIMOS Supernova Survey (SVISS). The sample 
contains 16 supernovae in total. The discovered supernovae have been classified as core collapse or type la supernovae (9 and 7, 
respectively) based on their light curves, colour evolution and host galaxy photometric redshift. The rates we find for the core collapse 
x 10~ 4 yr" 1 Mpc~ 3 h^ (with statistical and systematic errors respectively) at average redshift 0.39 and 
Mpc~ 3 hj at average redshift 0.73. For the type la supernovae we find a rate of 1 .29^57 x 10~ 4 yr -1 Mpc~ 3 
hL at (z) = 0.62. All of these rate estimates have been corrected for host galaxy extinction, using a method that includes supernovae 
missed in infrared bright galaxies at high redshift. We use Monte Carlo simulations to make a thorough study of the systematic effects 
from assumptions made when calculating the rates and find that the most important errors come from misclassification, the assumed 
mix of faint and bright supernova types and uncertainties in the extinction correction. We compare our rates to other observations 
and to the predicted rates for core collapse and type la supernovae based on the star formation history and different models of the 
delay time distribution. Overall, our measurements, when taking the effects of extinction into account, agree quite well with the 
predictions and earlier results. Our results highlight the importance of understanding the role of systematic effects, and dust extinction 
in particular, when trying to estimate the rates of supernovae at moderate to high redshift. 

Key words, supernovae - general, galaxies - distances and redshifts, galaxies - stellar content, surveys 



1. Introduction 

The cosmic rate of supernovae is an important observable that 
can be used to constrain the properties of galaxies at high red- 
shifts and to study the supernovae (SNe) themselves. Depending 
on the origin of the SN explosion, thermonuclear or core col- 
lapse, different aspects can be studied. T he first measure ment 
of the cosmic SN rate (SNR) was done bv lZwickvl dl938l) who 
found that "the average frequency of occurrence of supernovae 
is about one supernova per extra-galactic nebula per six hundred 
years" for the local volume. It is not until the latest decades that 
the higher redshift regimes have been possible to reach. 

More recently, large programmes targeting type la super- 
novae at intermediate and high (> 0.1) redshifts have been con- 
ducted to measure the expansion of the universe and do precision 
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cosmology (e.g. J Perlmutter et a lj 119971 IS chmidt et al. 1998; 
lAstier et all l2006t iRiess et alJ l2007t iMiknaitis et all ]20d7n . 
Some of these s urveys also report superno va rates for la SNe 
(e.g. iNeill et alJ 120061 iDahlen et all 120081) out to z ~ 1.5 
and core collaps e supernovae (CC SNe) (IDahlen et al] f2004: 
iBazin et alJl2009l) out to z ~ 0.7. Large surveys targeting SNe 
of any type have also been successful in finding and char- 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile, under ESO 
programme ID 167.D-0492. 



acterising supernovae as well as determining both local and 
intermediate redshift cosmic supernova r ates d Botticella et al.1 
2008; iDilday et al.ll20 1 Obt | Li et alJl201lB.ISharon et al l (120101). 
Dilda v et al.l (l2010al)T and iBarbarv et alJ J2012al) survey galaxy 
clusters, where the SN la rates are likely to be enhanced, to find 
supernovae and have r eported cluster SN rates out to redshift 0.9. 
IBarbarv et al.l d2012bl) also report SN la rates from detections in 
the foreground and background of the targeted galaxy clusters 
out to z ~ 1.5. 

The standard observational strategy for finding SNe at high 
redshift is to perform survey observations on a given field 
and then monitor the same field over a long period of time. 
Supernovae are discovered by searching the image s for vari- 
able sources using image subtraction tools (such as Alard 2000) 
to minimise subtraction residuals. The cadence of observations 
during the survey period is normally chosen to sample signature 
features of the SN light curves and colour evolution at the target 
redshifts. In this way photometric typing of the SNe is possible 
and the light curves can be used to study the SN characteristics. 
For the la surveys with cosmology as the main goal, follow-up 
spectroscopy of SN la candidates is needed to get a secure deter- 
mination of the redshift, to improve the accuracy in the distance 
measurement, and confirmation of the type. When calculating 
supernova rates from this kind of survey, care has to be taken to 
avoid selection effects from the spectroscopic observations. 
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Supernova typing normally includes studying the spectra of 
the SNe close to their peak luminosity and identifying spec- 
tral lines, notably H, He and Sill lines, something which is 
observationally very expensive and in practice unfeasible at 
high redshift for fainter SN types. Another method to type 
SNe is to compare the observed light curves and colour evo- 
lution to pre-existing templates of different SN types, i.e. pho- 
tometric S N typing. These methods have been demonstrated to 
work (e.g. iKuznetsova & Connoilvl2007HPoznanski et al.l2 007 ; 
iRodnev & Tonrvll2009l: lKessler et al.ll2010r) using somewhat dif- 
ferent techniques. Better typing accuracy is achieved with prior 
information on the redshift th rough photometric or spectroscopic 
redsh ift of the host galaxy (iKessler et al.ll2010t iMelinder et al.l 
l201lh . 

Thermonuclear supernovae (or SN la's) are thought to be 
white dwarfs that explode when they accrete matt er and ap- 
proac h the Chandrasekhar limit (for a review, see Leibundgut 



2000). When taking the luminosity-stretch relation (Phillips 



1993) into account the peak luminosity of these SNe exhibit 



a very narrow spread and can thus be used to accurately mea- 
sure cosmological distances. The exact details of the explosions 
and of the progenitor system are not fully understood. For ex- 
ample, the time between formation of the progenitor system 
and the supernova explosion - the so called delay time - is 
unknown. This delay time depen ds on the natu re of the com- 
panion star to the white dwarf (iGreggiol 12005). By studying 
the rates of la supernovae and comparing to either the cosmic 
star formation history (e. g. Dahlen & Fransson 1999; Neill et al. 
2006; Strolger et al.ll2010l) . or the star formation rates and stel- 
lar m a sses of the host ga laxies dSullivan et al. 2006; Totan feTaTI 
2008; Maoz et al. 2011) it is possible to set constraints on the 
delay time and thereby also on the progenitor system. 

Core collapse supernova explosions (CC SNe) are the end- 
points of the live s of mass i ve sta r s, with masses between ~8Mq 
and ~ 50 M Q (iNomotol [1981 iTsuiimoto etal] 119971; ISmarttl 
120091) . Since massive stars are short-lived compared to the cos- 
mic time-scales the CC SNe trace active star formation. By aver- 
aging the CC SN rate over cosmic volume the rest frame rate 
of star formation in that volume can be studied. In this way 
an independent measure of the star formation history at high 
redsh i ft can be obtained (iDahlen et alJ l2004t [ Cappellar o et all 
120051: iBotticellaet 10120081; iBazin et al.l 120091) . More conven- 
tional methods of finding the cosmic star formation rates in- 
clude measurin g the rest-frame UV light from galaxies at a given 
redshift (e.g. iGiavalisco et al.l2004tlBouwens et al.ll2009l) . mea- 
suring the far-infrared (FI R) light to take star for mation hidden 
by dust into account (e.g. iLe Floc'h et al] 20051) and deriving 
the ra tes from Ha measurements (Shi m et al]|2009HHaves et al] 
The UV and Har-based methods have a drawback, that is 
also present for the SNR method, in that a correction for dust 
extinction needs to be applied. Methods based on using the FIR 
light to estimate the added star formation from re-radiated UV 
light make it possible t o correct the star formatio n history for 
dust extinction effects. Hop kins & Beacoml d2006l) presented a 
compilation of the star formation history from multiple sources 
for z ~ - 6, taking the obscured star formation into account, 
and fitted an analytical function to the data. 

The light from a supernova has to pass through its host 
galaxy before starting on the long trip to reach our telescopes. 
When travelling through the gas and dust inside its host the su- 
pernova light will be subject to varying degrees of extinction, de- 
pending on the dust content of the galaxy and the position of the 
SN with respect to the observer (e.g. a SN situated in an edge- 
on galaxy will suffer from higher extinction, on average, than 



one in a face-on galaxvl. lHatano et all (Tl 998) and lRiello & Patad 
(2005) present models of extinction for core collapse and ther- 
monuclear SNe in normal spiral galaxies. These models are built 
by using Monte Carlo simulations of supernova positions within 
a galaxy with given morphology and dust content. By using the 
extinction models it is p ossible to estimate the effect on the ob- 
served supernova rates (IDahlen et al.l l2004). It should be noted 
that this method is mainly applicable to normal galaxies with 
low to medium amounts of gas/dust, typical of galaxies in the 
local volume of the universe. 

As the redshift increases, more and more of the star forma - 
tion takes place inside dust y galaxies. ILe Floc'h et al. (2005); 
iMagnelli etal] (120091 1201 ll) find that the star formation from 
these infrared bright galaxies dominate the global star forma- 
tion at redshift 0.7 and higher. In these galaxies the SN ex- 
plosions can be completely obscured by the large amounts of 
dust within the nuclear regions. For low to moderate amounts 
of extinction this can be estimated and taken into account in the 
light curve analysis. But for host galaxies with high dust content 
(such as luminous and ultra-luminous infrared galaxies, LIRGs 
and ULIRGs) most of the SNe may not be detectable, even in 
the n ear-infrared (NIR) where the extinction is strongly reduced 
(e.g. [Mattila & Meikle1 l2001l) . When calculating the rates, the 
number of missing SNe due to the change in average extinction 
in sta r forming galaxies w ith redshift needs to be compensated 
for dMannucci et al.ll2007l) . The derivation of the de-bias factors 
that are needed to correct the rates for this effect is further com- 
plicated by the recent discovery that the population of U/LIRGs 
at lo w redshift is quite different from the ones at high redshif t 
(e.g.lDaddi et al.l2010tlEibaz et al.l201lllKartaltepe et al.l201 lb . 
iMattila et al.l d2012l) have recently estimated the fraction of SNe 
missed in such galaxies both locally and as a function of redshift 
making use of the most recent results on the nature of U/LIRGs 
at different redshifts. 

The Stockholm VIMOS Supernova Survey (SVISS) is a 
multi-band (R+I) imaging survey aimed at detecting supernovae 
at redshift ~0.5 and derive thermonuclear and core collapse su- 
pernova rates. The supernova survey data were obta i ned o ver a 
six month period with VIMOS/VLT. IMelinder etal] (120081) de- 
scribe the supernova search method along with extensive testing 
of the image subtraction, supernova detection and photometry. 
The discovery and t yping of 16 supernova e in one of the search 
fields is reported in Mel inder et al.1 d201 ll) . Here we present the 
supernova rates estimated from the survey data along with delay 
time distributions for the la SNe and star formation rates (SFR) 
calculated from the CC SNe rates. 

The first part of the paper contains a description of the obser- 
vations and supernova sample. In Sect.[3]we describe the method 
used to determine the supernova rates and in Sect.|4]the resulting 
supernova rates, delay time distribution and star formation rates 
are presented. In the final section of the paper, Sect. [5] we dis- 
cuss the results and compare them to the work of other authors. 
The Vega magnitude system and a standard ACDM cosmology 
with {Ho,CI m ,Qa} = {70,0.3,0.7} have been used throughout 
the paper. 

2. The data 

2.1. Observations 

The data u sed in this paper wer e obtained with the VIMOS 
instrument dLe Fevre et al.l I2003I) mounted on the ESO Very 
Large Telescope (UT3) at several epochs during 2003-2006. The 
VIMOS instrument has four CCDs, each 2kx2.4k pixels with a 
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pixel scale of 0.205 "/pxl, covering a total area of roughly 4x56 
sq. a rcmin. The observatio ns were obtained in the ELAIS-S1 
field dLa Franca et al.ll2004h . in five broad band filters (U, B, V, 
R and /) centred at a = 00 /, 32'"13 s , 6 = -44°36'27" (J2000). 
The supernova search filters were R and /. Observations in these 
filters were divided into seven search epochs and 2 additional 
reference epochs (one before and one after). The search epochs 
were separated by roughly one month. The observational pro- 
gramme did not include any spectroscopy of the detected SNe 
and the analysis presented here is based on the R + 1 photometric 
data exclusively. 

The individual frames in each epoch were reduced, includ- 
ing removal of fringes, registered to a common frame of refer- 
ence and stacked. Each epoch image was photometrically cali- 
brated using photometric standard stars observed during one of 
the nights and local standard stars selected in the field. Detailed 
measurements of the seeing in each of the frames were done by 
modelling the point spread functio n. For more details on t he ob- 
servations and data reduction, see M encfa Trincha nt et a i] d20"l2h 
and lMelinder et all (120111) . 



Table 1. Supernovae in the SVISS 



SVISS ID 


Type 


P(Type) 


z 


SVISS-SN43 


la 


0.886 


0.43 


SVISS-SN161 


la 


1.000 


0.50 


SVISS-SN115 


la 


1.000 


0.40 


SVISS-SN116 


la 


1.000 


0.55 


SVISS-SN309 


la 


1.000 


0.47 


SVISS-SN402 


la 


1.000 


0.22 


SVISS-SN135 


la 


0.950 


0.98 


SVISS-SN14 


CC 


0.623 


0.36 


SVISS-SN51 


CC 


1.000 


0.51 


SVISS-SN54 


CC 


0.812 


0.77 


SVISS-SN261 


CC 


0.734 


0.37 


SVISS-SN55 


CC 


0.995 


0.83 


SVISS-SN31 


CC 


0.999 


0.12 


SVISS-SN56 


CC 


0.930 


0.57 


SVISS-SN357 


CC 


1.000 


1.40 


SVISS-SN24 


CC 


0.643 


0.81 



Notes. P(Type) is the probability (Bayesian evidence) for the best fit- 
ting main type (TN/CC). 



2.2. The supernova sample 

The supernovae were detected by using a dedicated pipeline 
that was developed and thoroughly tested by our team. For a 
detailed description along with the results of the testing see 
iMelinder et al.l (P2008h . The pipeline includes the following steps, 
first the reference image and the search i mage are sc aled to a 
common seeing using the ISIS 2.2 code {Alard 2000) and the 
reference image subtracted from t he search image. Auto matic 
source detection using SExtractor dBertin & Arnoutslll996h to- 
gether with a by-eye inspection is then done on the subtracted 
image to find variable objects. The final step is the aperture pho- 
tometry, and point spread function (PSF) modelling based pho- 
tometry, of the sources. Detailed simulations of the photometric 
accuracy are done to make sure that the error estimates are valid. 

To avoid including spurious detections in the output cata- 
logues we required that the supernova candidates were detected 
in both bands in two consecutive epochs. A detection is here de- 
fined as being brighter than the 3<x limiting magnitude in that 
epoch. Furthermore, we used the late control epoch to remove 
non-SNe (most likely active galactic nuclei, AGN) from the sam- 
ple, since no SNe are expected to be visible ~ 1 year after explo- 
sion at the wavelengths and redshifts considered here. It should 
be noted that the decision to require detection in both the R and / 
filters has a significant impact on the number of detected SN la's. 
These SNe are inherently redder than the CC SNe. At z > 0.8 
their light starts to become redshifted beyond the R filter and 
they are less likely to be detected in both of the filters used in 
this survey. This effect is taken into account in the rate calcula- 
tions (see Sect. [3}. 

For the sample of supernovae considered in this paper we 
used the pipeline described above on the ELAIS-S1 R + I data 
and found a total of 16 supernovae, seven la SNe and nine CC 
SNe, ranging from z ~ 0.1 to z ~ 1.4 (see Table [TJ. The light 
curves and colour evolution of the SNe were used to classify 
them into either la-like or CC types (9 different subtypes, see 
Table [2]). We then co-added the likelihoods for all subtypes be- 
longing to the type la or CC SNe, respectively, to find the most 
likely type. This was done using a Bayesian model selection 
code, where priors were used for the redshift (host galaxy photo- 
metric redshift), absolut e magnitude di stribution, extinction and 
time after explosion. In Meli nder et al.l d201 ll) we describe in de- 
tail how the sample was obtained together with the typing code. 



In that paper we also investigated the misclassification errors 
when applying our code to the SN sample by using a Monte 
Carlo simulated sample of ~ 18000 SNe and a local SN sample 
from SDSS-II containing 87 SNe. The resulting errors are used 
to estimate systematic errors for the SN rates (see Sect. 14.31 

Over the full redshift range where we have found SNe (z < 
1.4) we expect approximately 8% of the found la SNe, and 2% 
of the CC SNe, to be false positives. On the other hand, it then 
follows that only about 2% of the la SNe are mistyped compared 
to about 8% of the CC SNe. This systematic effect is thus not 
symmetric, overall we expect the la count to be slightly inflated 
and the CC count slightly deflated. 

2.3. Supernova redshifts 

The reference epoch R and / images were used together with 
the UBV deep images to obtain photometric redshifts for the 
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Fig. 1. The redshift accuracy for simulated supernovae with (up- 
per panel) and without (lower panel) a host galaxy photomet- 
ric redshift. The ~ 2000 SNe of all different subtypes used to 
make this figure are a sub-sample of the full sample of simu- 
lated light curves. The redshift on the x axis (z s ; m ) is the initial 
adopted redshift of the supernovae. The dotted lines indicate the 
Kfit £sim 1/(1 + W> > 0.15 limits. 
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supernova host galaxies. The redshi fts are calculated using the 
GOODSZ code (iDahlen et al.ll2010h . This is a x 2 template fit- 
ting code that uses em piri cal spectral energy d istributions from 
IColeman et all (Il980h and lKinnev et all (1 19961) . The photomet- 
ric redshifts were found to have an accuracy of cr- = 0.07 
where the redshift scatter Az is given by Az = cr z (l + z). 
D etails on the photometric redsh ift technique will be presented 
in Mencfa Trinchant et al.l (120121) . We require the galaxies to be 
detected in at least 3 filters to trust the resulting redshifts, some- 
thing which is fulfilled for all of the 16 supernovae. The hosts 
were selected by choosing the closest galaxy in physical dis- 
tance, calculating the distance from the photometric redshift and 
the angular distance. For two of the supernovae (SN309 and 
SN357) this means that the closest galaxy in terms of angular 
distance is not the host galaxy. These galaxies are at much higher 
redshifts (z > 3), thus it is extremely unlikely we would be able 
to detect a supernova in the m. More details on th e host galaxy 
identification is presented in lMehnder etaDd20Tlh . 

The photometric redshifts, together with the 68% confidence 
limits (as estimated by the x 1 fitting), were used to construct a 
Bayesian prior for the supernova typing. The typing code then 
outputs a most likely redshift for the best fitting supernova typ e 
for each supernova (details are given in iMelinder et al.ll201 ll) . 
This is the redshift estimate used for the supernovae in the rate 
calculations. Fig. Q] shows the redshift accuracy for the Monte 
Carlo simulated supernova sample described in Sect. 12.21 When 
typing the simulated supernovae we assume that the redshift er- 
rors for the host galaxies are equal to the redshift scatter, Az. The 
resulting redshift accuracy for the SNe with host galaxy redshifts 
is given by cr z = 0.07 with an outlier fraction (defined as having 
ksim - Zfltl > 0.15) of 3% over the full range of redshifts. The 
resulting redshift accuracy from the "photometric redshift+SN 
typing scheme is thus unchanged from the pure photometric red- 
shift scheme. 

For reference, we also show the accuracy for SNe typed with 
a flat redshift prior (i.e. sources without any photometric red- 
shift information). This is worse, with <j z =0.11 and an outlier 
fraction of 4%. For both populations it should be noted that the 
majority of the outliers comes from the z s im = 0.2 - 0.4 region 
and is due to mistyping of the supernovae. 

The typing was also rerun using a flat redshift prior for all 
of the discovered SNe. The fitted redshifts from this run were 
all within the 68% confidence limit of the host photometric red- 
shifts and none of the SNe changed main type. Based on this the 
chosen hosts are thus likely to be the correct ones, although it 
is certainly possible that some hosts have been misidentified. In 
any case, neither the typing nor redshifts are appreciably affected 
by the possible host misidentification. 

3. Supernova rate determination 

The sup ernova rates were determined using a Monte Carlo 
method (Dahlen et al. 2004). Using a set of supernova templates 
(see Table |2]l we simulate a number of supernova light curves 
of different types and at different redshifts with time of explo- 
sion set within our detection window. The requirement that each 
SN has to be detected in two consecutive epochs means that this 
window stretches from approximately 5 months before until five 
months after the start of monitoring. The main input parameters 
are the intrinsic supernova rates (la and CC) for a given redshift 
bin and volume-limited fractions of the different types. The tem- 
plate light curves ar e calculated from absolute magnitude light 
curves (mainly fr om Richard son et al.ll2002t IDahlen et al.ll2 004; 
iRichardson et al.ll2006l) and a set of spectra from lNugenj (120071) . 



Table 2. Properties of the supernova photometric templates 



Type 


M B 


0~M 


Stretch 


Fraction 


Ia-91T 


-19.64 


0.30 


1.04 


0.10 


la - normal 


-19.34 


0.50 


1.00 


0.52 


la - faint 


-18.96 


0.50 


0.80 


0.18 


Ia-91bg 


-17.84 


0.50 


0.49 


0.20 


Ibc - bright 


-19.34 


0.46 


N/A 


0.016 


Ibc - normal 


-17.03 


0.49 


N/A 


0.324 


IIL 


-17.23 


0.38 


N/A 


0.061 


Iln 


-18.82 


0.92 


N/A 


0.051 


IIP 


-16.66 


1.12 


N/A 


0.547 



Notes. (M B ) is the absolute magnitude in the lohnson-B filter at peak, 
o~m is the dispersion in the peak magnitude. Fraction refers to the frac- 
tion within each main type in a volume limited sample. 



For more information on how the template light curves are built 
see lMelinder et al.l d201 lb and references therein. 

The volume tric fractions for the la subtypes are adapted from 
Li et al.l d201 lb ), treating their Ia-HV sub-category as la-normal 
SNe. The la-faint category is used in the typing of supernovae 
to be able to better characterise normal la supernovae with a 
somewhat lower value of stretch. The fraction of low-stretch la 
SNe is estimated from the distribution of stre tch value in the 
Supernova Legacy Survey (ISullivan et al.ll2006l) la sample. The 
core collapse fract i ons ar e ba sed on a compil ation of supernovae 
from lSmartt et al.l (120091) and lLi et al.l(l2011bl) . treating lib SNe 
as Ibc-normal supernovae, since, with our time sampling of the 
light curves, these subtypes will look very similar. 

Each supernova is also given a host galaxy extinction. For 
type l a SNe we use the paramet risation of the Riello & Patatl 
(2005) simulations presented in Neil l et all d2006h . while for 
core collapse SNe we use Monte C arlo simulations based on 
the models of Riello & Patat (2005). The effects of spiral arms 
and the bulge component are considered negligible in the core 
collapse SN simulations and are not included. The extinction is 
scaled with the V-band optical depth through a simulated face- 
on galaxy at zero radius with Ty(0) = 2.5, which provides a 
reasonable match to the observed host galaxy extinct ion distri- 
bution for CC SNe within 12 Mpc dMattila et alJl2012l) . To cal- 
culate the wavelength dependence of the e xtinction, we use a 
Cardelli extinction law dCardelli et al.ll 19891) with Rv = 3.1 for 
Type la SNe and a Calzetti attenuation law dCalzetti et all 2000) 
with Ry = 4.05 for core collapse SNe. In Sect. 14.31 we inves- 
tigate how the choice of extinction models and extinction laws 
affects the results. It is important to note that this extinction is the 
result of simulations of normal spiral/elliptical galaxies. In this 
case the modelled mean extinction only goes above Ay ~ 1 mag. 
for galaxies with inclination higher than 60 degrees. According 
to the simulations, no supernovae in these normal galaxies ex- 
hibit very high extinctions. But observati ons of local supern ovae 
show that this is untrue, e.g. SN 200 2hh dPozzo et alJl2006l) and 
SN 2009hd dElias-Rosa et alJl201 ll) both have host galaxy ex- 
tinctions of Ay ~ 4 mag. The extinction adopted for our simu- 
lated supernova light curves thus only include the effect of small 
to moderate levels of extinction. In Sect. l3.2l we describe how the 
output rates are de-biased to account for the missing population 
of highly extinguished supernovae in normal galaxies as well as 
in U/LIRGs. 

The light curves are then fed through the same detection pro- 
cedure as the real supernova light curves (see Sect. 12.21 ). It should 
be noted that even with a quite conservative photometric limit 
(3cr in four data points) there are issues with completeness. The 
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detection efficiency of the survey does affect the observed rates 
since it starts to drop below one before the 3cr-limit is reached. 
The detection efficiencies described below are thus used to give 
a probability that a supernova is detected in a certain epoch and 
filter, and are included in the Monte Carlo simulations. We have 
used the detection efficiencies derived for hosts of intermediate 
brightness for all supernovae. The effect of this choice on the 
rates is small (see Sect. 14.31 ). 

The output number of detectable SNe from the simulations 
is then compared to the observed number in the chosen redshift 
bin. Since the redshifts for the SNe have quite high uncertainties 
(being based on photometric redshifts) we distribute the SNe in 
the bins according to a Gaussian distribution with the redshift 
uncertainty Az = 0.07(1 + z). We denote the number of redis- 
tributed SNe in each bin N re dist- The simulations are then iterated 
until the output number matches N re di s t- At that point the input 
intrinsic SN rate that produces the correct number of observed 
SNe is chosen as the true supernova rate. 

3. 1 . Detection efficiencies 

The detection efficiency for each epoch and filter was determined 
by simulating supernovae at variable brightness in the actual 
search images and then running the standard detection pipeline 
on the simulated fra mes. The procedures used are described in 
Meli nder et al.l d2008l) . We placed supernovae in different envi- 
ronments, i.e. with different host galaxy brightness, to study the 
effect of background light on the efficiencies. The detection ef- 
ficiencies are notably worse when the SNe are situated in bright 
galaxies (defined as having mi < 22). In Fig. [2] a sample of the 
detection efficiency testing is presented, showing the efficiencies 
for the three different host galaxy brightness modes. It should be 
noted that the 50% efficiency limits roughly corresponds to the 
3cr rejection limits. 

To use the efficiency curves in the Monte Carlo simulations, 
and to smooth out possible outliers in the measured efficiency, 
we parametrise the efficien cies using a S-curve parametrisation 
(as previously used by, e.g. lStrolger et al.ll2004l) given by: 



where e(m) is the fitted detection efficiency and m the magni- 
tude. £q is the maximum efficiency (which is very close to one 
for most of our fits), m c the magnitude when the efficiency has 
dropped to 50 %, and S a parameter that determines how fast the 
drop occurs. The best-fitting parameters for each epoch, filter 
and host galaxy brightness are found with a simple x 1 optimisa- 
tion algorithm and are shown as lines in Fig. [2] 

3.2. De-biasing the rates for extreme host galaxy extinction 

The star formation in the local universe takes place mostly in 
galaxies with low amounts of dust and thus the supernovae dis- 
covered in the local universe often have low extinctions. The 
use of core collapse supernovae as trace rs of recent star forma - 
tion is thus feasible at low redshifts (e.g. lBotticella et aTll2012l) . 
However, as we go to higher redshifts the bulk of the star forma- 
tion takes place in dusty galaxies with high infra red luminosities 
(LIRG s and ULIRGs ). A n umber of studies (e.g. lLe Floc'h et alj 
2005; Magnelli et al. 2009l l2011h have found the fraction of star 
formation taking place in LIRGs and ULIRGs to increase rapidly 
towards z ~ 1 , such that approximately half of the star formation 
at z ~ 1 is taking place in these galaxies. Furthermore, highly 




25.0 25^5 26^0 26^5 27.0 
R 




23.5 24.0 24.5 25.0 25.5 26.0 26.5 



I 

Fig. 2. An example of the detection efficiency for a given epoch 
for different host galaxy brightness. The efficiencies from simu- 
lations of supernovae in the field are given by the points while 
the lines are the S-curve fits to the efficiencies as described in 
the text. Errors for the measured efficiencies are binomial. Note 
that a total of 30 (one for each of five epochs x, two filters x, 
and three host galaxy types) S-curve fits have been used to esti- 
mate the detection efficiencies of the full survey. This figure is 
available in colour in the electronic version of the article. 

extinguis hed supernovae are present also in normal spiral galax- 
ies (e.g. Mattila et al1 l2012h . From now on we denote the effect 
these two factors have on the supernova rates obscuration, to dis- 
tinguish it from the standard host galaxy extinction of light from 
SN e in normal spiral and elliptical galaxies. 

iMannucci et alJ (l2007h compiled the star formation densities 
for different redshifts derived from UV and infrared observations 
in a number of studies. They used these results together with 
their own estimates on how many SNe are lost due to obscu- 
ration by dust in local starburst galaxies, LIRGs and ULIRGs 
to derive a correction for supernova rates at high redshifts. The 
estimates were bas ed on SN searches conducted in such galax - 
ies by that time dMaiolino et alj 120021: IMannucci et all 12003). 
which had found few SNe. Another caveat with the previous 
studies of this effect is that it was, at that time, almost com- 
pletely unknown what kind of galaxies the high redshift LIRGs 
and ULIRGs were. They were selected based on high luminosi- 
ties in t he mid- and far-infrared wavelengths, and the assumption 
made in lMannucci et al.l d2007l) was that they were the same kind 
of starbursting, irregular and compactly star forming galaxies as 
in the local universe. At the time this was not an unreasonable 
assumption, but later developments (Spitzer and Herschel obser- 
vations in particular) have shown that the high redshift U/LIRG 
population is dominated by disk g alaxies with more constant, al- 
though high, star formation r ates dDaddi et al]|2010l ; lEIbaz et all 
l201lHKartaltepe et alJl20lTh . the so called main sequence (MS) 
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galaxies. These galaxies have a higher content of gas and dust 
than local disk galaxies, but they do not exhibit the kind of com- 
pact star formation found in local U/LIRGs. A smaller fraction 
of the hi gh redshift U/LIRGs are more similar to the local coun- 
terparts dElbaz et al.ll201 UlKartaltepe et al.ll201 lb . 



3.2.1. De-bias factors for the core collapse supernovae 

For SNe in normal galaxies we have already taken into account 
the effect of low to moderate amounts of extinction, but, as 
previously mentioned, there are supe rnovae with e xtrem e ex- 
tinction also in normal disk galaxies. Mat tila et al.l d2012l) find 
that 15^jq% of the supernovae in these galaxies have extinc- 
tions significantly highe r than predicted by simulations follow- 
ing the recipe in iRiello & Patatl d2005l) . These are SNe that 
would likely be missed in dusty regions of the normal galax- 
ies by magnitude limited surveys. In the local universe almost 
all U/LIRGs are also chara cterized as starburst galaxies (e.g. 
iPereira-Santaella et alj |2010). In these galaxies less than ~ 20% 
of the CC SNe can be detected in optical searches and in some 
extreme cases, such as Arp 220, the entire SN population is ob- 
scured by th e large amounts of dust and can on ly be detected 
in radio (e.g. lParra et al.ll2007l) . lMattila et ail (120121) investigated 
the SN population of the nearby LIRG Arp 299 and found that 
83 " 



+9 



of SNe exploding in this galaxy will remain undetected 
by optical SN searches. They then assume that Arp 299 is rep- 
resentative of compact starbursting U/LIRGs and adopt this fig- 
ure as the missing fraction of supernovae in these galaxies. The 
missing fraction of SNe in high redshift U/LIRGs is highly un- 
certain. For the main sequence (i.e. non-starburst) galaxies little 
is known about the extinction, although they have been found 
to be disk-like and no t as compact as their local counterparts 
dKartaltepe et alj|201 lh . IMattila et all d2012l) assume a missing 
fraction of 37^ % for these galaxies (see Sect. l4.3l for a discus- 
sion on the systematic errors resulting from this decision). 

IMattila et al] d2012h calculate the de-bias factor as a func- 
tion of redshift by using the relative contributions to the cosmic 
star formation density of normal galaxies (defi ned as galaxies 
with L lR < 10 n L G ), LIRGs, and ULIRGs from iMaenelli etalJ 
d201lh and the relative contribu tions of the com pact and main 
sequence channels from Kartaltepe et al] d201 lh . Furthermore, 
they assume that Arp 299 only represents the U/LIRGs that are 
starbursting, i.e. lie more than three times above the specific star 
formation rate (sSFR) MS locus at high redshift. 

We adopt the de-biasing factors for th e CC SN rates in th e 
two redshift bins under consideration from lMattila et al.l d2012h . 
The final de-bias factors to correct the CC SN rates for obscura- 
tion in our two redshift bins is 1.34 at 0.1 < z < 0.5 and 1.52 
at 0.5 < z < 0.9. In Sect. l4.3l we give the errors on these factors 
and discuss the systematic errors on the rates resulting from the 
calculation. 



3.2.2. De-bias factors for the type la supernovae 

The previous section dealt with de-biasing the CC SN rate, 
but it should be noted that also la supernovae may be missed 
in SN surveys due to high extinction in star forming galaxies. 
Multiple studies have shown that the rate of la SNe is corre- 
lated with the current star formation rate in the host galaxies 
dMannucci et al.ll2006t ISullivan et al.l l2006). This indicates that 
there is a so called prompt cha nnel of thermonuclear explosion 
dScannapieco & Bif dsten 2005), or that the delay time distribu- 



tion e xtends to time intervals as short as -420 Myrs ( Maoz "et~aT] 
1201 lh . 

No la SN has ever been detected in the compact star form- 
ing cores of starburst galaxies. But, given that the extinction is 
extreme in such environments and that la SNe are faint in ra- 
dio, this may be the result of selection effects. However, the de- 
lay time distribution of la SNe has so far only been dete rmined 
down to a time difference of 420 Myrs (Maoz et al. 201 1), while 
the starburst pha se in U/LIRG starburs ts normally has a duration 
of ~ 100 Myrs (Marcil iacet al.ll2006l e.g.). After this time the 
starburst normally enters the post-starburst phase wher e the lack 
of gas makes the star formation stop (e.g. lGotol l2007). Most of 
the dense dust and gas shroud in a starburst core will likely have 
been disrupted and blown off in the post-starburst phase (e.g. 
iTremonti et al. 1 12007b . and any la SNe exploding after a delay 
time of > 100 Myrs will thus not experience higher extinction 
than a supernova in a normal star forming galaxy. 

Similar to what has been found for CC SNe, also a small 
fraction of the la SNe will be subject of extreme extinction in 
normal star forming disk galaxies. This fraction will depend on 
the delay time distribution and is likely significantly lower than 
the fractions found for CC SNe. At higher redshifts a fraction of 
type la SNe will explode in the main sequence channel of the 
U/LIRG population. These SNe will, as the CC SNe, suffer from 
higher extinction than the normal extinction models predict, but 
no estimates on this effect can be found in the literature and a 
detailed study is outside the scope of this paper. 

The f raction of la SNe misse d at high redshift has been esti- 
mated by Man nucci et al.l (l2007h . but suffers from the same prob- 
lems as the core collapse estimate in addition to systematic errors 
from the assumptions on delay time distribution and starburst 
lifetime (few details on the calculations are given in that paper) 
as described above. Because of the added uncertainties inherent 
in estimating a missing fraction for type la SNe we do not try 
to compute any missing fraction for the la SNe. Furthermore, it 
should be noted that any assumptions on the delay time distribu- 
tion may introduce a circular argument. To limit the circularity, 
and to avoid using a de-bias factor for which the systematic error 
is unknown, we assume that no la SNe will be missed because 
of obscuration (but note that we do apply a correction for the 
normal host galaxy extinction). In Section |4~31 we estimate the 
systematic error resulting from this assumption. 

4. Results 

4.1. Core collapse supernova rates 

The core collapse supernova rate in two redshift bins (0.1-0.5, 
0.5-0.9) is shown in Fig.[3]and in Table[3] The two bins do not 
include all of the SNe in the sample. By binning differently the 
total number can be increased, but this comes at the price of 
higher statistical uncertainty (with more bins) and higher sys- 
tematic uncertainties (with wider bins). Both the rates corrected 
for extinction/obscuration and the raw rates are shown in the ta- 
ble but the figure show the corrected values. The two sets of 
error bars in Fig.[3]show the statistical errors, and statistical and 
systematic errors added in quadrature, unless otherwise noted. 

In Table |4] and Fig. [3] we show a comparison of our rates 
with rates reported by ot her authors. At low redshift we plot the 
rate estimates of LOSS dLi et al.ll2011al) . which is the largest SN 
survey to date in the local volume with a total number of core 
collapse S Ne of 440. We a l so plo t the low redshift rates deter- 
mined bv ICappellaro et al] (Il999h . iBotticella et"ai] (120121) , and 
Mattila et al. (2012). At slightly higher redshifts we plot the re- 
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Table 3. S VISS supernova rates 



la Supernovae 


z 




0.3 < Z < 0.8 


(z) 




0.62 


N raw 




5 


Nredist 




4.92 


Ri a x 10~ 4 (no ext.) 




i ns+0.69 +0.23 
1A,o -0.49 -0.24 


Rk x 10- 4 (ext.) 




1 9Q+0.88+0.27 
• -0.57-0.28 


CC Supernovae 


z 


0.1 <z<0.5 


0.5 < z < 0.9 


(z) 


0.39 


0.73 




3 


5 


N re dist 


3.28 


3.27 


R cc x 10" 4 (no ext.) 
R cc x 10- 4 (ext.) 


l OQ + 1.72+0.97 
l.OJ _ ogg Q79 

9 4 6 +2.30+1.75 
1.33 -1.42 


9 cq+2.14 + 1.17 
-1.29 -0.80 

A 91+3.49+2.94 
-2.05-1.98 


R cc x 10~ 4 (ext.+obsc.) 


"1 9Q+3.08 + 1.98 
-1.78-L45 


f. 4f)+5.30+3.65 
-3.12-2.11 



Notes, (z) is the volume averaged redshift over the given redshift range. The supernova rates Ri a /cc are in units of yr~' MpcT 3 bi„ and given with 
and without corrections for extinction and obscuration. N raw is the raw number of supernovae per bin and N redist is the number when taking the 
redshift uncertainty into account. 



suits from lBotticella et ail d2008|). wh ich also includes the data 
presented in ICappellaro et alJ d2005l) . The second survey with 
rates based on a large number o f CC SNe (117 SNe) is the 
Supernova Legacy Survey (SNLS, Bazin et al. 2009) which pro- 
vides a data point with small error at z — 0.3. None of these 
surveys include de-biasing for obscuration. 

At highe r redshifts we plot t he rates from the GOODs de- 
termined by Dahle net alJ d2004h . which includes the effect of 
extinction in normal spiral galaxies, but not de-biasing for ob- 
scurat ion and the recen t measurements from the Subaru Deep 
Field dGraur et alJ201 ll) . which includes de-biasing for obscura- 
tion and extinction in normal galaxies. For comparison, the fig- 
ure also shows esti mates of the star fo rmat ion history from two 
differe nt sources, Strolger et al. (2004) and Hop kins & Beacon] 
(I2006J) . scaled to supernova rates (see Sect. 14.41 ). Both are cor- 
rected for extinction. 

The rate at z = 0.39 is consistent with the rates reported 
by other authors, even though different approaches to extinction 
correction is used. Our CC rate at high reds hift agrees well wit h 
the observations of Dahl en et alJ d2004l) and lGraur et alJ d201ll) . 
For a discussion on how our supernova rate measurements com- 
pares to the star formation history, see Sect. [5] 

4.2. la supernova rates 

The la supernova rate in one redshift bin (0.3-0.8) is shown in 
Fig. [4] and in Table [3] The rates in the table are given with and 
without extinction correction (no de-biasing to account for ob- 
scuration have been done on the la rate). The choice of using one 
bin for the la SNe is motivated by the low number of supernovae 
found at both low (z < 0.3) and high (z > 0.8) redshifts. By using 
the redshift bin given above we minimise the statistical errors. It 
should also be noted that the SNe outside the bin have not been 
excluded from the rate calculation. They contribute to the N re dist 
of the bin, because they have redshift probability distributions 
that stretches into the redshift interval. 

In Tableland Fig. [4] we show a comparison of our rate with 
rates reported by other authors. At low z we show the la rates 



from the LOSS dLi et al.l2011al) and ICappellaro et al.1 d!999l) . At 
slightly higher z we com pare with the rates ca lculated from the 
SDSS supernova search (Dildav et al. 2010b) a nd with the re- 
sults fr om the STRESS (Bott icella et al.ll2008l) iRodnev & Tonrvl 
(l2010h presented revised la rates from the IfA Deep Survey out 
to z = 1 .05 using new techniques, lower than those previously 
reported. The SNLS measured the la rate at z = 0.47 using a 
large sample of supernovae, at similar redshift lPain et al.l (f2002) 
presented some of the first rate s determined from c osmological 
survey data. At higher redshifts Dahle n et al.1 d2008) determined 
la rates from t he GO ODS (extending their sample from 2004). 
Barbarv et al. (2012b) has measured the volumetric SN la rate in 
the background gal axies of the Hubble Space Telescope Cluster 
Supernova Survey. Graur et aT| (1201 ll) reported the findings of 
the supernova search in the Subaru Deep Field, their rate mea- 
surements are the highest redshift measurements available and 
reach z ~ 1.7. 

Our rate measurement is higher than the results from many 
other surveys at this redshift, but marginally consistent with all 
other measurements when taking statistical and systematic er- 
rors into account. While the uncertainty of our rate measure- 
ment is too high and the redshift coverage too small to allow 
for detailed fitting of different DTD models, we can compare 
our rates to models introduced by other aut h ors dSullivan et all 
120061: iNeill et al.l 120061: fStrolger et all 120101: iGraur et al.ll201 ll) 
together with their chosen star formation history (see Sect. 14.51 
for details). 

4.3. Analysis of errors for the supernova rates 

The statistical errors are calculated using the prescription of 
iGehrelsl (|J986). The redshift bins have been chosen to provide 
a reasonable number of sources in each bin to get similar statis- 
tical errors in each bin. 

We now proceed to study the systematic errors of our rate 
estimates. Given that the total number of SNe is quite low, the 
statistical errors are high. One of the goals of this study is thus 
to find out whether any of the systematic effects can introduce 
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Table 4. Supernova rates in the literature 



<z> 


Rla/CC X 10- 4 


Reference 




IaS 


upernovae 


< 0.0066 (< 28 Mpc) 


> 0.35 


Smartt et al. (2009)* 


0.01 


24+ 07 

"•^ -0.07 


Cappellaro et al. (1999) 


< 0.014 (< 60 Mpc) 


f) 9^c+0.034f+0.043\ 
u.iUJ_ 0034 (._ 0043 ; 


Li etal. (2011) 


0.0375 




Dilday et al. (2010) 


0.1 


fl 9SQ+0.052/"+0.018\ 

"• iJ -0.044V-0.001; 


Dilday et al. (2010) 


0.15 


U,JU -0.034L0.005' 


Dilday et al. (2010) 


0.15 


A O9 +0.23/+0.07-) 
U,J -0.23^-0.06 ' 


Rodney & Tonry (2010) 


0.2 


f) -}4e+0.032/"+0.035\ 
u.jto_ 0030 v„ 0()7 ; 


Dilday et al. (2010) 


0.25 


3 65 +0 031 ( +0 182 ) 


Dilday et al. (2010) 


0.3 


n 454+0.037 (-+0.396X 


Dilday et al. (2010) 


0.3 


fl Q4+0.16/+0.2K 
KJ J -0.15^-0.22 > 


Botticella (2008) 


0.35 


n -34+0.19/+0.07A 
U - J -0.19^-0.03'' 


Rodney & Tonry (2010) 


0.442 


oo +0 - 50 r +00 °) 

-0.00 l -0.00*' 


Barbary et al. (2012) 


0.45 


3l +OA5 ( +OA2 ) 

U - J -0.15^-0.04 - 1 


Rodney & Tonry (2010) 


0.47 


fl 49+0.06^+0.13n 
U - H -0.06^-0.09 ' 


Neill et al. (2006) 


0.47 




Dahlen et al. (2008) 


0.55 


fl C4+0.099/+0.11 \ 
-0.086 '■-0.116' 


Pain et al. (2002)* 


0.55 


n Q9+0.14/+0.07-, 

"■^-o.i^-o.m' 


Rodney & Tonry (2010) 


0.62 


1 9Q+0.88/"+0.27\ 

-0.57M).28' 


SVISS (this work) 


0.65 


fl 4Q+0.17,-+0.14n 
^•^-0.17^-0.08' 


Rodney & Tonry (2010) 


0.74 


79+ - 33 


Graur et al. (2011), errors include systematics 


0.75 


fl /:o+0.21/+0.23\ 


Rodney & Tonry (2010) 


0.807 


1 1 Q+0.6O/+0.44-, 

-0.45^-0.28' 


Barbary et al. (2012) 


0.83 


1 3o +033 r +0 - 73 ) 

J -0.27 V- 0.51' 


Dahlen et al. (2008) 


0.85 


n 7 o+0.22c+0.3K 
^•/8_0.22Lo.l6-' 


Rodney & Tonry (2010) 


0.95 


76 +a25 ( +a32 ) 

-0.25^-0.26 ' 


Rodney & Tonry (2010) 


1.05 


n 7q+0.28,-+0.36n 
u - ' ^-0.28^-0.41 ' 


Rodney & Tonry (2010) 


1.187 


1 QQ+0.65C+0.69-, 
L - J:, -0A9^-0.26' 


Barbary et al. (2012) 


1.21 


1 Q9+0.36/+0.38-, 
1J -0.29^-0.32' 


Dahlen et al. (2008) 


1.23 


0.84^ 


Graur et al. (2011), errors include systematics 


1.535 


fl 77+I.O7 f+0.44\ 
U -' -0.5^-0.71' 


Barbary et al. (2012) 


1.61 


fl 49 +0.39/+0.19-, 
"•^-0.23 M).^ 


Dahlen et al. (2008) 


1.69 


1 - UZ -0.37 


Graur et al. (2011), errors include systematics 


CC Supernovae 


< 0.003 (< 11 Mpc) 


> 1 6 +04 


Botticella et al. (2012)*, stat. errors only 


- 0.003 (6-15 Mpc) 


> 1 5 +0 - 4 


Mattila et al. (2012)*, stat. errors only 


< 0.0066 (< 28 Mpc) 


> 0.88 


Smartt et al. (2009)* 


0.01 


0.58!™ 


Cappellaro etal. (1999) 


< 0.014 (< 60 Mpc) 


fl £9+0.067 /"+0.17\ 
U '° -0.067 ( --0.15-' 


Li etal. (2011) 


0.21 


1 14S+°- 43 C+ - 42 
la -0.34'-0.36 


Botticella et al. (2008) 


0.3 


1 49+0.3C+0.32-, 

-0.3 l --0.24'' 


Bazin et al. (2009) 


0.3 


9 cj+0.88/+0.75\ 


Dahlen et al. (2004) 


0.39 


-j 9Q+3.08/-+1.98\ 


SVISS (this work) 


0.66 




Graur et al. (2011), errors include systematics 


0.7 


3 96+ L03 (-+ L92 ) 


Dahlen et al. (2004) 


0.73 


6.40 ;:v( :f) 


SVISS (this work) 



Notes, (z) is the average redshift for the redshift interval, as given in the respective paper. The supernova rates Ri a /cc are in units of yr~' Mpc~ 3 
h 3 . Errors are statistical (systematical) unless otherwise noted. The rates are corrected for extinction given in the original reference. In a few cases 
(references marked with *) no extinction correction have been done. The local CC rates (papers by Botticella and Mattila) have been measured in 
small volumes and are likely dominated by a local overdensity of star formation. 
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Fig. 3. Core collapse supernova rates determined from the SVISS SNe (black stars). Grey stars are displaced by +0.01 in z and 
show the SVISS rates without extinction correction. Also shown is a collection of measured rates from other authors. We plot star 
formation histories from two different sources, scaled to supernova rates through the use of Eq. [2] Error bars are statistical with total 
errors (statistical and systematic added in qua drature) as a transparen t/faded extra error bar for all surveys. Redshift bin sizes are not 
shown, but are given in Table[3] The rate from Capp ellaro et alj d 19991) has been rescaled from SNu to a volumetric rate by assu ming 
a local B-band lumi nosity density 2.0 x 10 8 h 7 L Mpc" 3 at redshift ~ 0. The local CCSN rates from lBotticella et all d2012h and 
Matt ila et alJ d2012l) have been measured in small volumes and are likely dominated by inhomogeneity. The star formation rate in 
the local volume is higher than the average over an extended volume and can therefore cause the SNR to be higher than expected. 



errors larger than these. Calculating the mean statistical relative 
error (mean of the upper and lower limit differences) for the two 
types and all redshift bins we arrive at the following: 74% for 
CC at low z; 67% for CC at high z; and 56% for the la bin. These 
can be compared to the relative systematic errors calculated be- 
low. The second goal of the study is to make an extensive list of 
systematic effects that can affect supernova rate measurements 
and to try and estimate them. These effects will be very simi- 
lar for all supernova rate measurements of similar type and our 
compilation can thus be of use in future surveys which may be 
limited by systematic rather than statistical errors. A thorough 
study of the systematic effects is also of great help when trying 
to find better observational strategies in future surveys. 

A compilation of the systematic errors can be found in 
Table El 



Misclassification errors The systemat ic errors due to mis clas- 
sification in our sample are given in Melin der et al.l d201 lb and 
are on the order of 5-10 %. However, these estimates are really 
only valid when the number of detected SNe is large enough. For 
the small number statistic estimates presented here the system- 
atic errors have to be estimated in a different way. The expected 



numbers of misclassified SNe in each type and redshift bin is 
below one in all cases. To get a conservative estimate of the mis- 
classification systematic error we thus vary the observed number 
of SNe by one (positive and negative) and recalculate the rates. 
This test gives an error of approximately 20%, but is dependent 
on the actual number of sources in the given bin (a low source 
count gives a higher error). For bins with more than two SNe the 
error is smaller than the statistical errors. 

The systematic misclassification errors are overestimated 
when computed with this method, on average one SNe in our 
sample will have been misclassified, not one per bin and type 
as the estimate indicates. But it is impossible to say which bin 
and type that is affected by misclassification, hence we give this 
error estimate. 



Redshift The redshifts of the supernovae have uncertainties on 
the order of 0.1 when using the host photometric redshifts as 
a prior for the typing. This leads to a redistribution in redshift 
which affects the rates. To study this effect we redistribute the 
simulated SNe according to a Gaussian distribution assuming 
<x = 0.07(1 + z) and recalculate the rates. This test is repeated 
10,000 times to find the spread in the rates. The systematic error 



9 



I 

U 
CL 



5** 



-"J 1 

I 

O 



12 3 4 



J. Melinder et al.: The rate of supernovae at redshift 0.1-1.0 

Lookback time [Gyr] 



0.5 



0.1 



0.05 



5 

- r- 



8 

- r- 



9 

- r- 



10 

— I — 



11 




■ 


Dahlen 2008 


< 


Cappellaro 1999 


A 


Dilday 2010 


• 


Li 2010 


▼ 


Botticella 2008 


• 


Neill 2006 


♦ 


Graur2011 


► 


Pain 2002 (••■) 




Rodney & Tonry 2010 


X 


Barbary 2012 


★ 


SVISS 




HB06 SFH x 0.001 




Graur2011 DTD 




Stralger 2010 DTD 




Sullivan 2006 DTD 



0.0 



0.3 



0.6 



0.9 



1.2 



1.5 



1.8 



2.1 



2.4 



Z 



Fig. 4. Type la supernova rates determined from the SVISS SNe (black star). The grey star is displaced by +0.01 in z and show the 
rate without extinction correction. Also shown is a collection of published rates from the literature. We plot SN la rates resulting 
from assumptions on the delay time distribution combined with the star formation history. Error bars are statistical with total errors 
(statistical and systematic added in quadratu re) as a transparent/fad ed extra error bar for all surveys. Redshift bin sizes are not 
shown, but are given in Table[3] The rate from Cappe llaro et al.l d 19991) has been rescaled from SNu to a volumetric rate by assuming 
a local B-band luminosity density 2.0 x 10 h-jQ Lj sun Mpc -3 at redshift ~ 0. 



varies between ~15 to ~30% for the different types and redshift 
bins, lower than the statistical errors in all cases. 



Detection efficiencies We determine the systematic errors from 
the detection efficiency assumptions using the faint and bright 
host galaxy detection efficiency estimates. The binomial errors 
of the efficiencies are smaller than the difference between the 
results for the different host types in almost all cases. The effect 
on the rates is on the order of +2%, significantly lower than the 
statistical errors. 



Photometric errors The SVISS photometric zero point calibra- 
tion is accurate to within ~ 5%. Since we apply an absolute 3<x 
detection limit on our sample (see Sect. \2.2\ the estimated rates 
will be affected by a slight shift in the photometric zero points. 
This effect will in reality be random over the two filters and seven 
epochs, to get a conservative estimate we vary the detection lim- 
its by 0. 1 magnitude. The resulting rates vary by ~ 5%. 



Template choices Our selection of subtype light curve tem- 
plates and their assumed fractions influence the rates as well. 



The choices and assumptions we have made are based on obser - 
vational results dRichardson et al.ll2002l l2006t iLi et alJl20TTbh . 
The fractions are subject to statistical uncertainties as well as 
systematic uncertainties (e.g., evolution with redshift, selective 
obscuration). There is also the possibility that very faint super- 
novae are under-repr esented in magnitude-limited SN surveys 
dHoriuchi et al.ll2011l) . 

For the CC SNe we study the effect of introducing an ad- 
ditional light curve template, a IIP template with a peak Mb 
of -14.39 (as for the faint core collapse SN 1987A) with frac- 
tion that changes the percentage of CC SNe with Mb > -15 
from 7% to 30%. The effect of this change on the rates is quite 
large, it makes the rates go up by approximately 30%. This is 
smaller than the statistical errors but could help to explain any 
possi ble mismatch with the star formation history (as identified 
bv lHoriuchi et al. 2011). Varying the fractions randomly have a 
very small effect on the rates since the changes tend to cancel 
out. 

For the la SNe we have made tests with different combi- 
nations of using or not using the peculiar (91T and 91bg-like) 
and non-standard (la-faint) templates. We get the largest posi- 
tive change in the rates when putting the fraction of 9 IT SNe to 
zero (treating them as normal SN la's), on the order of 1-2%. 
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The largest negative change is found when setting the fraction 
of faint la supernovae (91bg and the la-faint subtypes) to zero, 
treating them as la-normal SNe, which results in a change of 
10-15% in the rates. These errors should be considered to be the 
extreme limits, since there is plenty of evidence that these faint 
la SNe do exist. 

Peak magnitudes As described above the choice of template 
distribution also affects the distribution of peak magnitudes for 
the simulated SNe. However, the average peak magnitudes used 
in this work have associated uncertainties (note that this is not 
the same as the scatter). 

For the CC S Ne the template peak m agni tudes have uncer- 
taintie s given bv iRichardson et alT(l2002l) and iRichardson ^tai] 
(2006). As conservative estimate we allow the peak magnitudes 
of CC SNe to vary by the lcr errors given in these papers and re- 
compute the rates. The resulting rates change by approximately 
10%, which is smaller than the statistical errors. 

For the la SNe the peak magnitudes are much better con- 
strained (due to their use in precision cosmology). Assuming an 
uncertainty of 0.05 magnitudes for the la subtype and 0.1 for 
the 91bg/91T subtypes we calculate the impact on the la rates. 
The effect is small, on the order of 3%, much smaller than the 
statistical errors. 

Host galaxy extincti on The effect of host g alaxy extinctions are 
based on models by Riello & Patatl 12005) for type la and core 
collapse SNe. For the type la SNe these models were used as 
given, while a different scaling was used for the CC SNe (see 
Sect. [3]). Using instead a positive Gaussian distributi on of E(B - 
V) with <tE(B - V)=0.2 mag. as in lNeill et al] d2006l) for Type la 
only changes the rate by ~ 2%. For core collapse SNe, we have 
examined the e ffect o f applying an extinction model following 
iRiello & Patal d200l as given. This has marginal effect on the 
rate in the low redshift bin (1%) and may affect the high redshift 

bin by 10%. 

Using the estimates from lHatano et al.l (Il998l) change the re- 
sulting rates by a very small amount (< 1%). We also study the 
effect of using different extinction laws for the extinction. For the 
CC SNe we use a Cardelli//?y = 3.1 law instead of the Calzetti 
law and for the la we use a st eeper Cardelli law with Ry = 2.1 
(as advocated b v lGoobarl2008l) . The impact on the rates is small, 
on the order of 1%. 

Dust obscuration in the host galaxy Manila etld] (120121) esti- 
mate lower and upper limits on the de-bias factors for CC SNe 
from different assumptions on the amount of obscuration in the 
SN host galaxies. The systematic errors on the CC SN rate re- 
sulting from these limits are: -10%/+28% at z = 0.39, and 
-10%/+33% at z = 0.73. The upper and lower limits on the 
de-bias factors corresponds to the upper and lower limits on the 
rates. 

For the la SNe we assumed a nominal de-bias factor of 1.0 
(i.e. a missing fraction of zero). As an upper limit to this value 
we have chose n to use a de-bias facto r calculated from the for- 
mula given in Ma nnucci et al ] (120071) for the redshift in ques- 
tion. The resulting systematic error on the la rates is then +8% 
at z = 0.62. 

AGN contamination We used a control epoch with observations 
obtained one year after our search period had ended to check 



whether any of our supernova candidates showed variability over 
this longer period. Except for the possibility of very rare peculiar 
SNe none of the SN subtypes should be detected one year later. 
The list of variable sources originally contained 3 1 entries, about 
ten of which showed variability over a long time-scale and were 
excluded from the final list of SNe. This routine is enough to 
limit the amount of AGN contamination to very low levels. We 
do not give this error in Table|5]since we estimate it to be smaller 
than all other systematic effects. 

Cosmic variance The size of the SVISS field is fairly small 
(4x56 sq. arcmin.). Cosmic variance is therefore a possible un- 
certainty in the rate estimat es. We have estimated th e cosmic 
variance by using the work of lTrenti & Stiaveftil (2008 )Q, and as- 
suming that the variance of SNe follow that of the overall galaxy 
population. The relative cosmic variance for the type la SNe red- 
shift bin is then 9%. For the CC SNe low redshift bin we estimate 
a relative variance of 12%, and 10% for the high redshift bin. 

In Table [5] we also give the sum of all the systematic errors 
(added in quadrature). These are also the systematic error esti- 
mates given in Table [3] The co-added systematic errors are all 
smaller than the statistical errors. The major contributions to the 
systematic errors comes from misclassification, template choices 
redshift uncertainties, and the obscuration correction. 

4.4. Comparing the CC supernova rate to the cosmic star 
formation history 

The core collapse supernova rate can be compared to the cosmic 
star formation history (SFH), since the time-scale for CC pro- 
genitors to explode is much shorter than the cosmic time-scale. 
The supernova rate is derived from the star formation history by 
assuming an initial mass function (IMF, denoted by f below) and 
the mass range of the SN progenitors (M\ and M u for the lower 
and upper mass limits, respectively). The volumetric rate of CC 
SNe, Rcc in units of yr~' Mpc 3 is then: 

R cc (z) = kxp»(z), (2) 

where p*(z) is the star formation history in units of Mq yr~' 
Mpc~ 3 . Rqc is a volumetric rate, and thus scales as h , this 
needs to be addressed when comparing to star formation histo- 
ries which may have a different scaling. The constant k is the per- 
centage of stars that explode as SNe per unit mass and is given 
by: 

_ C&WdM 

For this work we choose to use a Salpeter IMF (Salp eterl[l955h 
with progenitor masses between 8Mq and 50Mq. The constant 
k is then equal to 0.0070. 

The choice to use a Salpeter IMF is based on that most SFH 
measurements and other supernova rate measurements are us- 
ing this IMF. Recently a number of authors have argued that the 
Salp eter IMF is not consi stent with recent observations (see, e.g. 
Hop kins & Bea com 2006), and claim that an IMF with a flatter 
shape is better. The choice to use different IMFs introduce a sys- 
tematic uncertainty to the comparison. However, the IMF depen- 
dence is partly cancelled out by the dependence on IMF of the 

1 web based calculator available at 

http : //casa . Colorado . edu/~ trenti/CosmicVariance .html 
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Table 5. Systematic errors 



Error 


CC Supernovae 


la Supernovae 




0.1 <z<0.5 


0.5 < z < 0.9 


0.3<z<0.8 


Misclassification 


25% 


18% 


16% 


Redshift 


32% 


16% 


9% 


Detection eff. 


2% 


3% 


1% 


Photometric 


4% 


9% 


4% 


Template choices 


+31% 


+32% 


-12/+2% 


Peak magnitudes 


7% 


15% 


2% 


Extinction 


+ 1% 


+ 10% 


2% 


Obscuration 


-10/+28% 


-10/+33% 


+8% 


Cosmic variance 


12% 


10% 


9% 


Total systematic 


-44/+60% 


-33/+57% 


-23/+24% 


Mean statistical 


74% 


67% 


56% 



Notes. Errors are two-sided unless the sign is given. The total errors have been computed by co-adding the individual errors in quadrature. 



star formation history (where the IMF is used to normalise the 
star formation measurements ). To investigate this we compute k 
for the flattest IMF given in Hopkins & Beacom (2006) which 
is a modified Salpeter form from Baldrv & Glazebrook (2003) 
with a high mass slope of -1.15 (compared to -1.35 for the nor- 
mal Salpeter IMF). While the k factor changes to 0.0141 (a factor 
of ~ 2) this is countered by th e change of the SFH nor malisation 
by a factor of 0.5, given by Hopkins & Beacom (2006]). The to- 
tal change of the SNR-SFR comparison is on the order of 2%, 
significantly less than the statistical errors. Incidentally, this also 
tells us that there is little hope in trying to use the comparisons 
between SNR and star formation rates to constrain a possible 
universal IMF unless much larger supernova samples can be ob- 
tained. 

The upper limit for the progenitor mass comes from 
iTsuiim oto et al] d 19971) and is essentially the limit at which mas- 
sive stars will produce black holes instead of neutron stars with 
an associated supernova explosion. Varying the upper limit be- 
tween M u = [30, 125] M Q has a small effect on the SNR- 
SFR comparison, less than 10%, although this could be slightly 
higher with a flatter IMF. The lower limit o f 8 Mp) is grounde d 
on observations of supernova progenitors dSmartt et alj 12009). 
The choice of lower limit strongly affects the rate scaling but is 
luckily well constrained by these observations, varying the lower 
limit between their 68% confidence limits (7-9.5 Mq) changes 
k by +20%, also this lower than the statistical errors for our su- 
pernova rates. 



4.4.1. Star formation history parametrisations 

We use star formation histories, p*(z), from two di fferent sources 
as input for Eq. [2] [Hopkins & Beacom (2006) compile star 
formation measurements from many sources an d correct them 
for du st extinction based on the work done by iLe Floc'h et al.l 
(2005). The observations are then fit by the s imple analytical 
parametrisation introduced by [Cole et all (l200ll) : 



P*(z) 



(a + bz) h 



70 



1 + {zlc) d 



(4) 



where a, b, c and d are fitting parameters and hvo is the Hubble 
constant to account for changes to the a ssumed cosmology. The 
best fit from Hopkins & Beacom (2006) gives parameter values: 
a = 0.0170, b = 0.13, c = 3.3 and d = 5.3. They give a con- 
servative error estimate by plotting the confidence regions of the 
fit, at z ~ 1 the uncertainty is on the order of 30%. 



IStrolger et al.l (|2004|) fit the data compilation of 

iGiavalisco et al.l (|2004) by using a parametrisation of the 
form: 

pjt) = a(t h e-" c + de d{, -'° )lc ), (5) 

where a, b, c and d are parameters, fis the age of the universe 
and fo the age of the universe at z = (13.47 Gyrs with our 
chosen cosmology). The best fit parameter values are: a = 0.182, 
b = 1.26, c = 1.865 and d = 0.071. To compare this with our 
results we convert it to p*(z). 

4.5. Delay time distributions for la supernovae 

The delay time distribution (DTD) for la supernovae is the sub- 
ject of controversy in current research. By convolving an as- 
sumed DTD with the underlying star formation history the la 
rate can be calculated. This can then be compared with the ob- 
served rate s to put constraints on the DTDs. The la rate can thus 
be written (IStrolger et al]l2004l) : 



flla(0 



V f p*{t')<&{t- 

Jl 



t')dt\ 



(6) 



where v is the number of SNe formed per unit stellar mass, p*(f) 
the star formation history and <5(r) the delay time distribution 
which represents the percentage of supernovae that go off at time 
r after a single burst of star formation. Also, to is the time at 
which stars start to form in the universe, we assume this hap- 
pens at z = 10; corresponding to a time of 0.45 Gyrs after the 
big bang (with our chosen cosmology). The integral can easily 
be converted to redshift space to yield R\ a (z) by performing a 
variable substitution. 

As noted above, the uncertainties in our rate estimates make 
it quite fruitless to try and fit different DTDs to our data. Instead 
we choose to compare our rates to the best-fitting DTDs from 
the literature. In Fig. |4] we show four different choices of DTD. 
The simplest one is just a rescaling of the star formation history 
assuming that one la SN explode per IOOOMq of stellar mass 
formed dNeill et aT]l2006l) . This corresponds to using a Dirac 
delta function DTD with the peak at r = 0, i.e. no delay time 
between the SN explosion and star formation or, in other words, 
a pro mpt channel fo r the la pro genitor-to-supernova process. 

In lStrolger et all d2010l) and Vstrolger et al] d2004l) a Gaussian 
shaped (or close to Gaussian) DTD with a mean delay time of 
~ 3 Gyrs is found to give the best fit to the GOODs la supernova 
rates. We use the unimodal ske w-normal DTD param etrisation 
and best fitting parameters from IStrolger et all (2010) and con- 
volve this with the star formation history given in Eq.|5](the same 
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SFH as the one used by these authors). The DTD parametrisation 
is: 

4>(t) = — e~' 12 dt', (7) 

WW J-iaf 

with the best fi tting parameters a> = 0.2, £ = 3.2 and a = 2.2. 

iGraur et all (1201 ll) fit the results from the Subaru Deep Field 
supernova search together with the rates of many other surveys 
using different power law DTDs. The DTD models and SFHs 
they try show quite similar fitting quality, the rate we show in 
Fig. |4]uses a power law DTD with an exponent (3 equal to 0.97, 
their best fitting value: 



<D( T ) = OiA 



(8) 



where Oi is a normalisation parameter. <I> is also set to for 
r < 40 million years. This is then convolved with the SFH given 

in Eq.|4]to get the rates plotted in th e Fig. [4] 

IScannapieco & Bildsten (2005), Mannucc i et al. I d20Q6h and 
Sulli van et al.l (120061) find strong evidence for the existence of 
a prompt channel for la SNe by looking a t the properties of 
the host galaxies of la SNe (and other data). Ma oz et al. 
study the hosts of the local SNe found in the LOSS and find a 
prompt channel (with r < 420 Myrs) with 99% significance. 
Models of la explosions also show the p ossibility of such a 
channel dNomoto et al.l 120091: iRuiter et all 1201 lb . The la rates 
are then modelled with one prompt component, directly propor- 
tional to the SFH, and one dela yed component . We p lot the re- 
sulting la rate from the work bv lSullivan et al.l (]2006), who use 
a parametrisation of the form 



X 

Jo 



RuV)=A pAt')dt' + BpJt). 



(9) 



Note that the SFH used by us is not strictly the same as the one 
used by these authors and can thus not be used direc t ly in t his 
equation. We plot the SNR as given in ISullivan et all (120061) in 

Fig. a 



5. Discussion and summary 

We have presented supernova rates from the SVISS along with 
a description of the methods used to compute them and an ex- 
tensive analysis of systematic errors for the rates. The result- 
ing rates for the core collapse SNe are: 3 .29 + ,™ + !;?8 x 10" 



T OQ+j.OS+l.io y 

J " i -1. 78-1.45 A 

yr _1 Mpc~ 3 hL (with statistical and systematic errors respec- 
tively) at <z) = 0.39 and 6.40^0+3.65 x 10 -4 yr -i Mpc -3 ^ 
at (z) = 0.73. The CC rates have been corr ected for obscuratio n 
in dusty environments using the results of iMattila et al.l (120121) . 
as described in Sect. | 3.2land host galaxy extinction based on 
the models in iRiello & Patatl (120051) . Uncorrected values can be 
found in Table [3] 

The rate e s timate s follow the star formation history well. 
iHoriuchi et al.l d201 lb point out that the core collapse supernova 
rates found in both local and high redshift searches seem to be 
too low by a factor of two when compared to th e star formation 
history compiled in iHopkins & Beacom (2006). Our results do 
not show this difference, both of our rate estimates are consistent 
with the star formation history within the statistical errors. But 
our errors are quite large, the rates estimated by other authors, in 
particular at low z, do in fact differ significantly from the SFH. 
The strongest e vidence for this comes from the rate estimate by 
iLi et al.l (1201 lal) at low z which has small errors and is clearly be- 
low the SFH. At higher redshift the problem is less severe, which 



could be due to the increased statistical errors. At high redshift 
the problem may also be somewhat alleviated by the obscura- 
tion corrections, which we have included in our plotted rates in 
Fig. [3] the rates from other surveys plotted in this figure do not 
include this corre ction (with the exception of the data point from 
IGraur etal.ll20TTh 



Horiuchie taLl(l201ll) suggest that taking missed SNe due to 



extinction and dust obscuration in LIRGs and ULIRGs into ac- 
count is not enough to explain the difference. Instead they sug- 
gest that the reason is that the assumed fractions of faint and 
very faint CC SNe are too low. In our tests of the systematic 
uncertainties we find that assuming 30% of the CC SNe to be 
faint (Mb > -15) boosts the rates by ~30%, not enough to 
bridge the factor of two gap found for the local SN searches. 
Of course, the assumption on fractions of faint CC SNe may be 
very different for the different surveys, making it possible that 
other data points may go up more than this. While the SFHs we 
compare with have also been corrected for dust and hidden star 
formation, the supernovae are probably affected differently by 
the presence of obscuring dust. The number of supernovae that 
are missed in dusty starburst galaxies and LIRGs is currently not 
well constrained even in the local universe (e.g Ma nnucci et al.l 
120031: IMattila et 31.112 0041 The de-bias factors for extinction and 
dust obscura t ion in normal galaxies and U/LIRGs derived by 
Matti la et al.l (120121) can make the difference between the pre- 
dicted and observed CC SN rates disappear. In this paper we 
have used these factors to de-bias the core collapse supernova 
rates. The de-bias factors from this study are slightly larger 
at low redshifts, but lower at high redshift, than the factors 
from Mannucci et a II (12001 . The uncertainties - both statistical 
and sy stematic - of the de-bias factors derived in Mattila et al. 
(2012) have been thoroughly studied, and will hopefully de- 
crease as more observations of SNe in dusty galaxies are ob- 
tained. High angu l ar resolution observations at near-IR (e.g. 
Matti la et all l2007t iKankare et alJl2008l 120121) and rad io (e.g, 
iPerez-Torres et alj|2010t iRomero-Cafiizales et al. l l20Tll) wave- 
lengths have recently been used to detect and characterise the 
hidden SN populations in the nearest LIRGs. Such studies are 
needed to constrain the complete rates, properties and extinction 
distributions towards the CC SNe buried in such dusty galaxies. 
Eventually, these studies will hopefully provide a robust estimate 
for the numbers of SNe missed by optical searches both locally 
and at high-z. 



The resulting rate for the la supernovae is: 1.29 



+0.88 +0.27 
-0.57 -0.28 



10 yr Mpc h^ at (z) — 0.62. This rate has been corrected 
for host galaxy extinction, but not for high redshift obscuration 
(a positive correction for this is included in the systematic error). 

Because of the quite high statistical errors and the lack of la 
SNe beyond z — 1 we do not try to fit any DTDs to our la rate 
measurements. The comparison to the plotted models show that 
all these models are consistent with our rates. The high rate of 
la SNe at z ~ 0. 5 measured by SVIS S is in strong agreement 
with the results of Dahlen et al. (2008), and thus in support of a 
Gaussian-like, fairly wide, DTD for la SNe. However, it should 
be noted that our measurement is also consistent with the distri- 
butions proposed by other authors. The type la SN rate measure- 
ments at z > 0.4 differ by up to a factor of ~ 2. We believe that 
the cause of this large scatter is the large statistical and system- 
atic errors that the SN surveys suffer from in this redshift range. 
It should also be noted that there are differences in the methods 
used to calculate the rates in the different surveys. It is therefore 
of utmost impo rtance t h at the systematic errors are correctly es- 
timated. Strol ger et alj d2010l) find that models with a prompt 
la component are hard to reconcile with the rates measured at 
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redshifts higher than one. Increased sample siz es at these high 
redshifts, or more studies of la host galaxies dGallagher et alJ 
120051: ISullivan et alj|2006h are needed to constrain the contribu- 
tion from this channel. 

In Section l3.2.2l we outline our motivation for not de -biasing 
the type la SN rates to account for extremely high extinction 
in U/LIRGs. It is very important to understand that the missing 
fractions of type la SNe situated in U/LIRGs similar to Arp 299 
very likely depend on the DTDs. It is likely that also the missing 
fractions in other types of galaxies will depend on the DTD, al- 
though to a lesser extent. Using the type la SN rates to estimate 
the DTD thus inherently suffers from circularity to some extent. 
Assuming a zero missing fraction is a choice in itself - the DTD 
implicit to such an assumption has a cut-off at some delay time 
t < 200 Myrs. If the DTD has no such cut-off, and a significant 
number of type la SNe have very short delay times, our assump- 
tion is faulty. In that case the resulting rates may be off by more 
than the adopted systematic error. 

The determination of supernova rates at high redshift is dif- 
ficult. The SNe detected at high redshift will only be sampling 
the bright end of the SN luminosity function. This means that 
any global statistics estimated from such measurements will be 
sensitive to assumptions on the luminosity function made during 
the calculations. There are a number of additional assumptions 
going into the rate calculations that affect the rates differently. It 
is important to estimate the systematic errors these assumptions 
give rise to. With the exception of the misclassification error, 
that essentially scales with the sample size for samples smaller 
than ten given a misclassification ratio of 10%, we have shown 
that the systematic errors are on the order of 50% when using 
photometric redshifts and with the present uncertainties in tem- 
plate fractions and peak magnitudes. Furthermore, the assump- 
tions made when correcting the rates for extinction/obscuration 
are shown to have a large effect on the systematics. Presently 
little is known about the number of SNe missed in LIRGs and 
ULIRGs, especially important at high redshift, which is evident 
in the large uncertainties on the de-bias factors. 

Given the low numbers of SNe for most high redshift sur- 
veys it is perhaps tempting to try and use the rates found in dif- 
ferent surveys together when comparing to models (for the la 
SNe) and other sources (SFH for CC SNe). However, it's not 
straight-forward to do this. Different surveys estimate system- 
atic errors in different ways, some include more sources and 
some less. If combined as given using co-added statistical and 
systematic errors, the risk is that greater weight is given to sur- 
veys where fewer systematic error estimates are included (given 
that the sample sizes are similar). We believe that the work pre- 
sented in this paper shows the importance of including a variety 
of systematic effects to correctly estimate the uncertainties of su- 
pernova rates at high redshift. This will be even more important 
for future surveys with larger sample sizes and therefore lower 
statistical errors. 
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